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ABSTRACT 

In this paper, we extend previous work on the absorption of high energy 
7-rays in intergalactic space by calculating the absorption of 10 to 500 GeV 
7-rays at high redshifts. This calculation requires the determination of the 
high-redshift evolution of the intergalactic starlight photon field, including its 
spectral energy distribution out to frequencies beyond the Lyman limit. To 
estimate this evolution, we have followed a recent analysis of Fall, Chariot & Pei, 
which reproduces the redshift dependence of the starlight background emissivity 
obtained by the Canada-France redshift survey group. We also include the UV 
background from quasars. 

We give our results for the 7-ray opacity as a function of redshift out 
to a redshift of 3. We also give predicted 7-ray spectra for selected blazars 
and extend our calculations of the extragalactic 7-ray background from 
blazars to an energy of 500 GeV with absorption effects included. Our results 
indicate that the extragalactic 7-ray background spectrum from blazars should 
steepen significantly above 20 GeV, owing to extragalactic absorption. Future 
observations of a such a steepening would thus provide a test of the blazar origin 
hypothesis for the 7-ray background radiation. We also note that our absorption 
calculations can be used to place limits on the redshifts of 7-ray bursts; for 
example, our calculated opacities indicate that the 17 Feb. 1994 burst observed 
by EGRET must have originated at z <~ 2. Finally, our estimates of the 
high-energy 7-ray background spectrum are used to determine the observability 
of multi-GeV 7-ray lines from the annihilation of supersymmetric dark-matter 
particles in the galactic halo. 
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1. Introduction 

The EGRET experiment aboard the Compton Gamma Ray Observatory has detected 
more than 50 blazars extending out to redshifts greater than 2 (Thompson, et al. 1996). 
It is expected that 7-rays from blazars with energies above the threshold energy for 
electron-positron pair production through interactions with low energy intergalactic photons 
will be annihilated, cutting off the high energy end of blazar spectra. Such absorption 
is strongly dependent on the redshift of the source (Stecker, De Jagcr & Salamon 1992). 
Stecker & De Jager (1997) have calculated the absorption of extragalactic 7-rays above 0.3 
TeV at redshifts up to 0.54 and presented a comparison with the spectral data for the low 
redshift blazar Mrk 421. 

The study of blazar spectra at energies below 0.3 TeV is a more complex and physically 
interesting subject. In addition to intergalactic absorption, one must be able to distinguish 
and to separate out the effects of intrinsic absorption and natural cutoff energies in blazar 
emission spectra. Initial estimates of intergalactic absorption of 10 to 300 GeV 7-rays in 
blazar spectra at higher redshifts have been given by Stecker (1996), Stecker & de Jager 
(1996) and Madau & Phinney (1996). However, in order to calculate such high-rcdshift 
absorption properly, it is necessary to determine the spectral distribution of the intergalactic 
low energy photon background radiation as a function of redshift as realistically as possible. 
This calculation, in turn, requires observationally based information on the evolution of the 
spectral energy distributions (SEDs) of IR through UV starlight from galaxies, particularly 
at high redshifts. Conversely, observations of high-energy cutoffs in the 7-ray spectra 
of blazars as a function of redshift, which may enable one to separate out intergalactic 
absorption from redshift-independent cutoff effects, could add to our knowledge of galaxy 
formation and early galaxy evolution. In this regard, it should be noted that the study of 
blazar spectra in the 10 to 300 GeV range is one of the primary goals of a next generation 
space-based 7-ray telescope, GLAST (Gamma-ray Large Area Space Telescope) (Bloom 
1996) as well as a number of ground-based 7-ray telescopes currently under construction. 
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2. Redshift Dependence of the Intergalactic Low Energy SED 

Our main goal is to calculate the opacity of intergalactic space to high energy 7-rays 
as a function of redshift. This depends upon the number density of soft target photons 
(IR to UV) as a function of redshift, whose production is dominated by stellar emission. 
To evaluate the SED of the IR-UV intergalactic radiation field we must integrate the total 
stellar emissivity over time. This requires an estimate of the dependence of stellar emissivity 
on redshift. Previous work has either assumed that all of the background was in place at 
high redshifts, corresponding to a burst of star formation at the initial redshift (Stecker 
1996; Stecker & De Jager 1996; MacMinn and Primack 1996) or strong evolution (similar to 
a burst), or no evolution (Madau and Phinney 1996). In this paper, we use a more realistic 
model which is consistent with recent observational data. 



2.1. Basic Calculation of Stellar Emissivity 

Pei & Fall (1995) have devised a method for calculating stellar emissivity which bypasses 
the uncertainties associated with estimates of poorly defined luminosity distributions of 
evolving galaxies. The core idea of their approach is to relate the star formation rate 
directly to the evolution of the neutral gas density in damped Lya systems, and then to use 
stellar population synthesis models to estimate the mean co-moving stellar emissivity S,^{z) 
(erg/s-cm^-Hz) of the universe as a function of frequency u and redshift z (Fall, Chariot & 
Pei 1996). Our calculation of stellar emissivity closely follows this elegant analysis, with 
minor modifications as described below. 

Damped Lya systems are high-redshift clouds of gas whose neutral hydrogen surface 
density is large enough (> 2 x 10^° cm~^) to generate saturated Lya absorption lines in the 
spectra of background quasars that happen to lie along and behind common lines of sight 
to these clouds. These gas systems are believed to be either precursors to galaxies or young 
galaxies themselves, since their neutral hydrogen (HI) surface densities are comparable to 
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those of spiral galaxies today, and their co-moving number densities are consistent with 
those of present-day galaxies (Wolfe 1986; see also Peebles 1993). It is in these systems that 
initial star formation presumably took place, so there is a relationship between the mass 
content of stars and of gas in these clouds; if there is no infall or outflow of gas in these 
systems, the systems are "closed" , so that the formation of stars must be accompanied 
by a reduction in the neutral gas content. Such a variation in the HI surface densities of 
Lya systems with redshift is seen, and is used by Pei & Fall (1995) to estimate the mean 
cosmological rate of star formation back to redshifts as large as z — 5. 

Pei & Fall (1995) have estimated the neutral (HI plus Hel) co-moving gas density 
Pc^g{z) in damped Lya systems from observations of the redshift evolution of these systems 
by Lanzetta, Wolfe, & Turnshek (1995). (Here Pc = Sifg/STrG is the critical mass density 
of the universe. The deceleration parameter is assumed throughout to be Qq = 0.5, with 
cosmological constant A = 0.) Lanzetta, et al. have observed that while the number density 
of damped Lya systems appears to be relatively constant over redshift, the fraction of higher 
density absorption systems within this class of objects decreases steadily with decreasing 
redshift. They attribute this to a reduction in gas density with time, roughly of the form 
0,g{z) — Qpoc^, where Pc^go is the current gas density in galaxies. Pei & Fall (1995) have 
taken account of self- biasing effects to obtain a corrected value of fig{z); we have reproduced 
their calculations to obtain flg{z) under the assumptions that the asymptotic, high redshift 
value of the neutral gas mass density is Qg^i — 1.6 x 10~^/io^, where ho = ifo/(100 
km/s-Mpc). In a "closed galaxy" model, the change in co-moving stellar mass density 
PcCls{z) = —pcClg{z), since the gas mass density Pc^g{z) is being converted into stars. This 
determines the star formation rate and consequent stellar emissivity (Pei & Fall 1995). 

To determine the mean stellar emissivity from the star formation rate, an initial 
mass function (IMF) 0(M) must be assumed for the distribution of stellar masses M in a 
freshly synthesized stellar population. To further specify the luminosities of these stars as a 
function of mass M and age T, Fall, Chariot, & Pei (1996) use the Bruzual-Charlot (BC) 
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population synthesis models for the spectral evolution of stellar populations (Bruzual & 
Chariot 1993, Chariot & Bruzual 1991). In these population synthesis models, the specific 
luminosity Lstari^^, M,T) (erg/s-Hz), of a star of mass M and age T is integrated over a 
specified IMF to obtain a total specific luminosity Si,{T) per unit mass (erg/s-Hz-g) for 
an entire population, in which all stellar members are produced simultaneously {T — 0). 
Following Fall, Chariot, and Pei (1996), we have used in our calculations the BC model 
corresponding to a Salpeter IMF, 0(M) dM oc M'^-^^ dM, where O.IMq < M < 125Mq. 

The mean co-moving emissivity S,y{t) is then obtained by convolving over time t the 
specific luminosity with the mean co- moving mass rate of star formation, pcflf- 

S,{t) = pc tdt'S,{T = t - 0^/(0- (1) 
Jo 

Note that the star mass formation rate pcflf{t) that appears in this equation is not the 
same as PcCls{t), the change in total stellar mass density. This is because Pc^s is the 
rate at which mass is permanently being converted into stars; since some stellar mass is 
continuously being returned to the interstellar medium (ISM), the instantaneous mass rate 
of star formation Pc^f is larger than Pc^ls, the two being related by 

Clsit) = Qfit) - f dt' R{t - t')nf{t'), (2) 
Jo 

where R{T), provided by the BC models, is the fraction of the initial mass of a generation 
of stars formed at T = that has been returned to the ISM. 

2.2. Metallicity Corrections 

The BC models' specific luminosities S^{T) are calculated assuming that the metaUicity 
content Z during star formation is fixed at our current solar metallicity value {Zq = 0.0169). 
However, the metallicity content of the universe is not static, but evolves with redshift 
as early populations of stars return freshly synthesized metals to the interstellar medium 
during their various phases of mass loss. For example, in a survey of 1/3 of the known 
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damped Lyman-alpha absorbers, Pettini et al. (1994) found that the typical metaUicity is 
0.1 that of the present solar value at a redshift oi z ^ 2. Since the specific luminosity of 
a star of a given mass is also a function of its metaUicity content (lower metallicities give 
bluer spectra), the metaUicity of a stellar population must be taken into account when 
integrating the mean emissivity over redshift. 

The effect of metaUicity content in stellar population models has been examined 
by Worthey (1994). Using the IMF 0(M) oc M-'^-^^dM with O.IMq < M < 2Mq, 
Worthey has calculated the mass-to- light ratios {M / L) as a function of population age T 
and metaUicity Z, for the color bands U through M. We have plotted his {M/L) values for 
the U and B bands in Figures 1 and 2 respectively. One can see that for a fixed metaUicity, 
the logarithm of the luminosity decreases approximately linearly with the logarithm of 
population age, and that for a fixed age the U and B luminosities decrease as the metaUicity 
increases. We have made a hnear fit to each fixed- metaUicity {M/L) (T) computed data 
set, obtaining a metaUicity correction factor factor Cx{Z/Zq) — Lx{Z)/Lx{Zq), where X 
designates the color band. ^From the parallel linear fits made to the computed data for each 
y = log{Z/Zo) value, it is seen that a common correction factor, C, applied to each y ^ 
data set will bring these data into rough agreement with the y — values of {M/L) (T). 
These correction factors are plotted in the inset figures, whose abscissa is Z/Zq and whose 
ordinate is the correction factor. Our fit to Worthey's computed data in Figure 1 gives a 
continuous correction factor JC{Z) for A = 0.35 //m, the center of the U band. Similar fits 
to the B band data (Figure 2) and to the V, /, and K band data (not shown) result in the 
relation 



log£(A,Z) = log 



L{X,Z) 
L{X,Ze) 



0.33- 



0.30 



A 



log 



+ 



■"3y 



0.066 - 



0.063 



A 



log 



Z^ 



0. 



(3) 



for 0.3 /xm < A < 2.2 /im. Outside of this wavelength region we take 

£(A < 0.3/xm,Z) = jC{X = 0.3/xm,Z) and £(A > 2.2//m, Z) = £(A = 2.2//m,Z). 
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Note that increased metallicity gives a redder population spectrum (Bertelli et al. 1994). 

Limitations to this correction factor include the fact that Worthey's calculations only 
apply to stars with ages greater than T = 1.5 Gyr, and that the upper mass limit of his IMF 
{2Mq) is much lower than that of the BC model which we employ (125Mq). Additional 
uncertainty exists below 0.3 yum since Worthey's calculations extend only to the U band. 
We have chosen to assume a constant enhancement factor below A = 0.3/im. For all of 
the above reasons, our enhancement factor C is really a conservative lower limit to the 
corrections to the BC models in the ultraviolet. 

Population synthesis models in which varying metallicity is included do exist (Bertelli 
et al. 1994), and efforts to reconcile differences in computed spectra generated by these 
various models have been made (Chariot, Worthey, and Bressan, 1996). 

2.3. Absorption by Dust and Gas within Source 

The emissivity given in Eq. |l| assumes that all stellar emission escapes from the 
gas system which contains the stars. However, some absorption of stellar radiation occurs 
both by dust and gas within the larger damped Lya systems. Above the Lyman limit, 
this absorption is dominated by dust, while below the Lyman limit, absorption by neutral 
hydrogen and singly-ionized helium dominates. Defining the mean transmission fractions, 
averaged over the optical depths of damped Lya systems, by 7dust(i^, 2;) and T^as^v.z), the 
final expression for the effective stellar emissivity is 

—^-—-^nfiz')Ciu, z')S4T = t{z) - t{z')]. (4) 

The distribution of optical depths of Ly« clouds due to dust is can be adequately 
represented by f{x) = fox^^'e'^, where x = t^/t^^z), t^{z) being a characteristic (redshift 
dependent) cloud dust opacity, and a 1 (Fall et al. 1996). Under the assumption that 
both dust and stars are uniformly distributed throughout each Lya cloud, the fraction of 
radiation T^^\Td) produced by stars in a given cloud of optical depth that escapes dust 
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absorption is given by 




2 1 - e"^*? 

TT] [{l + r]) + (1 -ri)e-^'^\ ' 



(5) 



where r] = vT^^TJ, and u is the average albedo of dust, taken to be the same value as in 
our Galaxy {u ~ 0.4 to 0.6; Whittet 1992). (We have calculated Eq. | using the 2-stream 
approximation [Chandrasekhar 1950]). We note that the dust opacity in Eq.^is assumed 
to be proportional to the HI surface column density A^hi and metallicity Z, 



where -E(i^) is the normalized galactic interstellar dust extinction curve (Savage and Mathis 
1979). Integrating Eq.^ over the Lya opacity distribution function / of Pei & Fall (1995), 
we obtain 7dust('^, z), and find it to have a minor effect on the emissivity, 7dust('^, z) being 
typically of order unity. 

Below the Lyman limit (A < 0.0912 /im), the opacity is dominated by neutral gas 
absorption: Tgi^u) = A^hiO"hi('^) + ^HeiO"Hci(^), where cthi and UHei are the HI and Hel 
photoionization cross sections (Osterbrock 1989). With the A'^hi and A^Hei distributions of 
the Lya systems being related to the dust opacity distribution /(r^) through Eq.^, the 
distribution for Tg can be obtained. Integrating Eq.|| (now with r] = 1), weighted with the 
Tg distribution, gives Tgas{h',z). 



Figure 3 shows the calculated stellar emissivity as a function of redshift at 0.28 yum, 
0.44 /im, and 1.00 /im, both with and without the metallicity correction factor C{X, Z). We 
have also plotted the observations of the cosmic emissivity by the Canada-French Redshift 
Survey (Lilly, Le Fevre, Hammer, & Crampton 1996) at these rest-frame wavelengths for 
comparison. With a lower mass cutoff of O.IMq in the IMF, we obtain emissivities which 
are roughly a factor of 2 higher than those obtained by Lilly, et al. (1996). To bring our 
emissivities down to the observed values requires that we reduce the lower mass limit in 
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2.4. 



Numerical Results 



the IMF to O.O2M0, which puts a fraction (0.45) of the mass into effectively nonluminous 
compact objects. We note that a similar reduction was achieved by Fall, et al. (1996) by 
modifying the power law index in the IMF; a higher index results in a lower emissivity (Pei 
1996, personal communication). 

Overall, our emissivities, both with and without the metallicity corrections, are in 
reasonable agreement with the data at lower redshifts (Lilly, et al. 1996). Although the 
differences for £1, between the no-metallcity and metallicity cases for 2; < 1 are not great, 
they become substantial at larger redshifts for both optical and UV wavelengths. This has 
notable effects on the opacity of the radiation background to high energy 7-rays , as will be 
seen in Section 4. We note that our dotted-line curves in Figure 3 (no metaUicity correction) 
are essentially a reproduction of the emissivities calculated by Fall, et al. (1996). 

In all cases as shown in Figure 3, the stellar emissivity in the universe peaks at 
1 < z < 2, dropping off at both lower and higher redshifts. Indeed, Madau, et al. (1996) 
have used observational data from the Hubble Deep Field to show that metal production 
has a similar redshift distribution, such production being a direct measure of the star 
formation rate. (See also the review by Madau (1996).) 



The co-moving radiation energy density Uu{z) (erg/cm'^-Hz) is the time integral of the 
co-moving emissivity Siy{z), 



where u' = v{l -\- z')/{l -\- z) and z^^;^ is the redshift corresponding to initial galaxy 
formation. The extinction term e^'^"^ accounts for the absorption of ionizing photons by the 
clumpy intergalactic medium (IGM) that lies between the source and observer; although 
the IGM is effectively transparent to non-ionizing photons, the absorption of photons by 
HI, Hel and Hell can be considerable (Madau 1995). The presence of damped Lya and 



2.5. Calculation of the Diffuse Radiation Energy Density 
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Lyman-limit systems (Lanzetta, et al. 1995) and the Lyman-alpha forest, coupled with the 
absence of a HI Gunn-Peterson effect (Gunn & Peterson 1965; Steidel & Sargent 1987) 
indicates that essentially all of the HI, Hel, and Hell exists within intergalactic clouds 
whose measured HI column densities range from approximately 10^^ to 10^^ cm~^. 

The effective optical depth Tefj between a source at redshift z' and an observer at 
redshift z owing to Poisson-distributed intervening Lyman-alpha clouds is given by (Paresce, 
McKee, & Bowyer 1980) 

rM.,...')^l"d."iyN,,^(l-e<-y), (8) 

where r(z/) = [Nuicrm{i^) + iV'HoiO"Hoi(z^) + A^HeiicrHeii(z^)], i^" = ^^(1 + z)/{l + z"), and 
is the distribution function of clouds in redshift z and column density A^hi- As pointed 
out by Madau & Shull (1996), when r <^ 1, r^s is just the mean optical depth of the 
clouds; when r ^ 1, Tcfr becomes the number of optically thick clouds between the source 
and observer, so that the Poisson probability of encountering no thick clouds is e~'^<=* , as 
required. 

For the distribution function of Lyman-alpha clouds we use the parameterization of 
Madau (1995) (see also Miralda-Escude & Ostriker, 1990, Model A2): 

d^N { 2.4 X lO'^N^l^l + z)2-4^ 2 X 10^2 < jy^^^ < ^ 59 ^ lO^W^ 

(9) 



9Nmdz I ^ g ^ io8^r-/-5(i + ^)0.68 1,59 ^ iQi^ < Nui < lO^^cm-^. 

Using Eqs. ^ and ^ and the stellar emissivity S,^{z) in Eq.|^, we obtain the background 
energy density Uu{z), shown in Figures 4 and 5, calculated with and without the metallicity 
correction, C, respectively. These also include the contribution to the UV background from 
QSOs (Madau 1992), which are believed to dominate the diffuse background radiation 
below the Lyman hmit and to be responsible for the early (z > 5; see Schneider, Schmidt, 
& Gunn 1991) reionization of the IGM. 

Although it is possible that UV emission from QSOs alone may be able to account 
for the nearly complete reionization of the IGM (Meiksin & Madau 1993; Fall & Pei 1993; 
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Madau & Meiksin 1994), it has been argued that additional sources of of ionizing radiation 
are required (Miralda-Escude & Ostriker 1990), these perhaps being young galaxies which 
leak a fraction (up to ~ 15%) of their ionizing radiation through HII "chimneys" (Dove & 
ShuU 1994; Madau & ShuU 1996). We have therefore assumed in our calculations that 15% 
of the stellar emission escapes from the galaxies (protogalaxies) through these chimneys, 
unattentuated by dust or gas. (We note, however, that recent observations of four starburst 
galaxies by the Hopkins UV Telescope (Leitherer et al. 1995) indicate that less than 3% of 
Lyman continuum photons escape from these sources.) Figures 4 and 5 indicate that in our 
calculation the A < 0.0912/xm background is indeed dominated by QSOs, so that the actual 
value of the escape fraction we choose is not too significant. 

The intergalactic energy densities given in Figures 4 and 5 are quite consistent with 
the present upper limits in the UV (Martin & Bowyer 1989; Mattila 1990; Bowyer 1991; 
Vogel, Weymann, Ranch & Hamilton 1995). It should be noted that our results as shown 
in figures 4 and 5 give emissivities from starlight only and do not include dust emissivities 
in the mid-infrared and far-infrared. 

3. Opacity of the Radiation Background and its Effect on Blazar Spectra 

With the co-moving energy density u^{z) evaluated, the optical depth for 7-rays owing 
to electron-positron pair production interactions with photons of the stellar radiation 
background can be determined from the expression (Stecker, et al. 1992) 



where Eq is the observed 7-ray energy, Ze is the redshift of the 7-ray source, x = (1 — cos^^), 
being the angle between the 7-ray and the soft background photon, h is Planck's constant, 
and the pair production cross section cr-y-y is zero for center-of-mass energy ^/s < 2meC^, 
being the electron mass. Above this threshold. 



t{Eq, Ze) i dz— dx- di'{l + z) 



Jo dz Jo 2 Jo 



/■^e , dt , X 



3 

hv 



a^^[s = 2Eohux{l + zf], (10) 




(11) 



- 13- 

where /3 = (1 - ^mlc" j sfl"^ . 

Figures 6 and 7 show the opacity t(£^o, ^) for the energy range 10 to 500 GeV, 
calculated with and without the metallicity correction. Extinction of 7-rays is negligible 
below 10 GeV. Above 500 GeV, interactions with photons with wavelengths of tens of \vca. 
become important, so that one must include interactions from infrared photons generated 
by dust reradiation (Stecker & De Jager 1997), which we have neglected here. For 300 GeV 
7-rays , at redshifts below 0.5, our opacities agree with the with the opacities obtained by 
Stecker & De Jager (1997). Our calculated opacity, even with the metallicity correction, is 
probably somewhat low in the 10 to 30 GeV energy range, because we have underestimated 
the value of L in the UV (see previous discussion) . 

Note that these calculated opacities are independent of the value chosen for /iO) ^ seen 
in Eqs. 4, 7, and 10. The emissivity Ei, in Eq. 4 scales as /iq, since neither S^,^ C nor dtQ,f 
depends on h^, while pc scales as h^. Eq. 7 shows then that Ui, scales as h^, and in Eq. 10 
this Hq factor is cancelled by the integration over time t. 

With the 7-ray opacity t{Eq,z) calculated out to 2; = 3, the cutoffs in blazar 7-ray 
spectra caused by extragalactic pair production interactions with stellar photons can be 
predicted. Figure 8 shows the effect of the intergalactic radiation background on a few of 
the 7-ray blazars ("grazars") observed by EGRET, viz., 1633+382, 3C279, 3C273, and Mrk 
421. We have assumed that the mean spectral indices obtained for these sources by EGRET 
extrapolate out to higher energies attenuated only by intergalactic absorption. Observed 
cutoffs in grazar spectra may be intrinsic cutoffs in 7-ray production in the source, or may 
be caused by intrinsic 7-ray absorption within the source itself. Whether cutoffs in grazar 
spectra are primarily caused by intergalactic absorption can be determined by observing 
whether the grazar cutoff energies have the type of redshift dependence predicted here. 

Figure 8 indicates that the next generation of satellite and ground-based 7-ray 
detectors, both of which will be designed to explore the energy range between 10 and 300 
GeV, will be able to reveal information about low-energy radiation produced by galaxies at 
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various redshifts and at different stages in their evolution. 

4. Constraints on Gamma-ray Bursts 

Our opacity calculations have implications for the determination of the origin of 7-ray 
bursts, if such bursts arc cosmological. As indicated in Figure 6, 7-rays above an energy of 
~ 15 GeV will be attenuated if they at emitted at a redshift of ~ 3. On 17 February 1994, 
the EGRET telescope observed a 7-ray burst which contained a photon of energy ~ 20 
GeV (Hurley, et aL1994). If one adopts the opacity results which include our conservative 
metallicity correction (Figure 6), this burst would be constrained to have originated at a 
redshift less than ~2. (An estimated redshift constraint of ~ 1.5 was given by Stecker and 
De Jager (1996), based on a simpler model.) Future detectors may be able to place redshift 
constraints on bursts observed at higher energies. 

5. The High Energy Gamma Ray Background from Blazars 

In a previous paper (Stecker & Salamon 1996), we presented a model for calculating 
the extragalactic 7-ray background (EGRB) due to unresolved grazars. We gave results for 
7-ray energies up to 10 GeV (where there is effectively no 7-ray absorption) which were 
compared to preliminary EGRET data (Kniffen et al. 1996) Using the intergalactic 7-ray 
opacities calculated here, we can now extend the results of this EGRB model out to an 
energy of 0.5 TeV. 

Our EGRB model assumes that the grazar luminosity function is related to that of 
flat spectrum radio quasars (FSRQ), so that we can use FSRQ luminosity and redshift 
distributions (Dunlop and Peacock, 1990) to obtain a grazar luminosity function. The 
effects of grazar flaring states, 7-ray spectral index variation, and redshift dependence are 
also been included in this model; see Stecker and Salamon (1996) for details. By integrating 
the grazar luminosity function weighted by our new opacity results, we obtain a grazar 
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background spectrum up to 500 GeV which properly includes the effect of 7-ray absorption. 

Figure 9 shows this EGRB spectrum compared with the preliminary EGRET data. 
Note that the spectrum is concave at energies below 10 GeV, reflecting the dominance of 
hard-spectrum grazars at high energies and softer-spectrum grazars at low energies; it then 
steepens above 20 GeV, owing to extragalactic absorption by pair-production interactions 
with radiation from external galaxies, particularly at high redshifts. Both the concavity 
and the steepening are signatures of a blazar dominated 7-ray background spectrum. 

Because the extragalactic 7-ray background in our model is made up of a superposition 
of lower-luminosity, unresolved grazars, its intensity is determined by the number of sources 
in the Universe which are below the detection threshhold of a particular telescope. A 
telescope with a superior point source sensitivity gives a higher source count, thereby 
reducing the number of unresolved sources which constitute the diffuse 7-ray background. 
In Figure 9, the upper spectra which are close to the EGRET data are obtained using the 
EGRET threshold; the lower curves correspond to the projected sensitivity of the proposed 
next generation GLAST sateUite detector, which is expected to have a detection threshhold 
of ~ 2 X lO-'' cm-^s-i above 0.1 GeV. 

It should also be noted that above 10 GeV, blazars may have natural cutoffs in their 
source spectra (Stecker, De Jager & Salamon 1996) and intrinsic absorption may also 
be important in some sources (Protheroe & Biermann 1996). Thus, above 10 GeV our 
calculated background flux from unresolved blazars, shown in Figure 9, may actually be an 
upper limit. 

6. Observability Above Background of a MuIti-GeV Gamma-ray Line 

Produced by Neutralinos 

The nature of the dark matter in the universe is one of the most important fundamental 
problems in astrophysics and cosmology. The non-baryonic mixed dark matter model 
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with a total Q — 1 (Shafi & Stecker 1984) gave predictions for fluctuations in the cosmic 
background radiation (Schaefer, Shafi & Stecker 1989; Holtzman 1989) which were found to 
be in good agreement with the later COBE measurements. The best agreement appears 
to be found for ~ 20% hot dark matter, of which massive neutrinos are the most hkely 
candidates, and ~ 80% cold dark matter (Pogosian & Starobinsky 1993, 1995; Ma & 
Bertchinger 1994; Klypin, et al. 1995; Primack, et al. 1995; Liddle, et al. 1996; Babu, 
Schaefer & Shafi 1996). 

The most popular cold dark matter particle candidates are the lightest sypcrsymmetric 
particles (LSPs), the neutralinos (hereafter designated as x particles). Cosmologically 
important x particles must annihilate with a weak cross section, {(tv)a ~ 10~^^ cm^s~^; 
calculations show that such cross sections lead to a value for ~ 1 with VLy,o:{av)^ . 
The fact that supersymmetry neutralinos are predicted to have such weak annihilation 
cross sections is an important reason why they have become such popular dark matter 
candidates. Prehminary LEP 2 results give a lower limit on the mass of the x of ^ 20 
GeV (Elhs, Falk, Ohve & Schmitt 1996). 

In the minimal supersymmetry model (MSSM), x can be generally described as a 
superposition of two gaugino states and two Higgsino states. Grand unified models with 
a universal gaugino mass generally favor states where x is almost a pure B-ino {B) [e.g. 
Diehl, et al. 1995), but other states such as photinos and Higgsinos are generally allowed 
by the theory. Kane & Wells (1996) have presented possible accelerator evidence from CDF 
that X i^ay be a Higgsino of mass ~ 40 GeV. 

Dark matter neutralinos will produce 7-rays by mutual pair annihilation. This process 
is expected to occur because neutralinos are Majorana fermions, i.e., they are their own 
antiparticles. Indeed, most of this mutual annihilation would have occured in the very early 
universe, a process which determines the present ( "freeze out" ) value of Vl^ and leads to the 
relation oc (o"^w)~^, where the bracketed quantity is the thermal-averaged annihilation 
cross section times velocity (Ellis, et al. 1984). This leads to the relation that the 7-ray flux 



from neutralino annihilation is inversely proportional to Q-^^. Thus, the annihilation 7-ray 
flux is limited from below by cosmological constraints on the maximum value of 

There are two types of 7-ray spectra produced by XX annihilations, viz., (1) 7-ray 
continuum spectra from the decay of secondary particles produced in the annihilation 
process, and (2) 7-ray lines, produced primarily from the process xX ~^ 77 Rudaz 
1989). The cosmic 7-ray flux from xX annihilation is proportional to the line-of-sight 
integral of the square of the x particle density times {av)A- 

The continuum 7-ray production spectra from xX annihilation can be calculated 
for different types of neutralinos by starting with the appropriate branching ratios for 
annihilation into fermion-antifermion pairs which produce hadronic cascades leading to 
the subsequent production and decay of neutral pions (Rudaz & Stecker 1988; Stecker 
1988; Stecker & lyika 1989). Stecker & Tylka (1989) discuss in detail the various channels 
involved in continuum 7-ray production via xX annihilation and give the resulting spectra 
for some lower mass x particles. Such continuum fluxes from xX annihilations would be 
difficult to observe above the extrapolated cosmic background which we show in Figure 
9. However, with good enough sensitivity and energy resolution, it might be possibile to 
observe a two-photon annihilation line from xX annihilation. The general considerations 
for observability of this line were discussed by Rudaz & Stecker (1991). We update this 
discussion here, using (1) our new calculation of the 7-ray background flux from blazars 
shown in Figure 9, (2) recent accelerator limits on supersymmetric particle masses, and 
(3) the proposed sensitivity and energy resolution of a next generation space based 7-ray 
telescope taken from the G LAS T propossl (Elliot 1996). 

The energy of the xX ~^ 77 decay line is = M^. The line width is given by Doppler 
broadening. For galactic halo particles, this width is roughly Px^x ~ ^^~^^x^ n^uch 
smaller than the energy resolution proposed for any future 7-ray telescope. Upper and lower 
limits on Q yield lower and upper limits on the 7-ray line flux respectively (see above). 
Other limits can be obtained in flux-energy space (Rudaz & Stecker 1991). Accelerator 



determined lower limits on My, give lower limits on the line energy. Lower limits on the 
mass of the sfermion exchanged in the annihilation process give upper limits on {av)^ 
since {av)A oc Mj'^. In fact, since the particle density n^- = p^-fMy. and {av)A oc M^/Mj, 
the predicted annnihilation line flux (f){Ej) oc Mj"^. Further limits are obtained from the 
inequahty > M^, which is the tautology following from the condition that x be the LSP. 

If we assume that annihilations occur mainly through slepton exchange, i.e., Mg ^ Mf, 
we can obtain an upper limit on the 2j line flux. This is because LEP 1.5 gives a lower 
limit of ~ 70 GeV on the slepton mass (de Boer, Miquel, Pohl & Watson 1996), whereas 
the substantially higher squark mass lower limit of ~ 150 GeV would imply much lower 
fluxes, since 0-y oc MJ"^. 

The lower limit on the slepton mass implies an upper limit on the line flux from BB 
annihilation such that the event rate for a next generation 7-ray telescope with an aperture 
of 1 m^sr would be about 5 photons per year for a line in the energy range between 20 and 
100 GeV. (If the x particles are Higgsinos, the event rate would be much lower.) Using 
all of these constraints, the allowed region for a neutralino annihilation line in flux-energy 
space is plotted in Figure 10. 

In constructing Figure 10, we have used the GLAST proposed estimate of the 
point source sensitivity after a one-year full sky survey to estimate the background from 
unresolved faint blazars (see Figure 9 and the discussion in the previous section). We then 
obtain the background photon number for an appropriate exposure factor of 1 m^yr-sr and 
energy resolution of 10%, and plot the square root of this number, which represents the 
natural background fluctuations above which a line must be observed. Of course, a higher 
exposure factor would reduce the point source background and increase the sensitivity to a 
line flux, as would a better energy resolution. It should also be noted that the background 
above 10 GeV shown in Figures 9 and 10 may be overestimated (see previous section). 

Another possible way in which dark matter may produce 7-rays and neutrinos is if 
the LSP is allowed to decay to non-supersymmetric, ordinary particles. Supersymmetry 
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theories involve a multiplicative quantum number called R-parity, which is defined so that it 
is even for ordinary particles and odd for their supersymmetric partners. Thus, if R-parity 
is conserved, as is usually assumed, the LSP is completely stable, making it a potential 
dark matter candidate. However, such may not be the case. R-parity may be very weakly 
broken, allowing the LSP to decay with branching ratios involving 7-rays and neutrinos 
{e.g., Berezinsky, Masiero & Valley 1991). For x particles to be the dark matter, their 
decay time should be considerably longer than the age of the universe. 

The possible radiative decay X ^ + 7 will give a 7-ray line with energy = 
Such a line has the potential of being more intense than the annihilation line. Whereas the 
XX annihilation rate and consequent line fiux is cosmologically limited by requiring fi^ to 
be a significant fraction near 1 (see previous discussion), the decay-line fiux is limited only 
by the particular physical supersymmetry model postulated and constraints from related 
accelerator and astrophysical data. Thus, invocation of x decay involves a higher order of 
particle theory model building and speculation. We only wish to mention here that there 
is a possibility that a decay line may be sufficiently intense to be observable above the 
background. 

7. Conclusions 

We have calculated the 7-ray opacity as a function of both energy and redshift for 
redshifts as high as 3 by taking account of the evolution of both the SED and emissivity 
of galaxies with redshift. In order to accomplish this, we have adopted the recent analysis 
of Fall, et al. (1996) and have also included the effects of metallicity evolution on galactic 
SEDs. We have then considered the effects of the 7-ray opacity of the universe on 
7-ray bursts, blazar spectra, and on the extragalactic 7-ray background from blazars. In 
particular, we find that the 17 Feb. 1994 EGRET burst probably originated at 2; < 2. 
Because the stellar emissivity peaks between a redshift of 1 and 2, the 7-ray opacity which 
we derive shows little increase at higher redshifts. This weak dependence indicates that 
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the opacity is not determined by the initial epoch of galaxy formation, contrary to the 
speculation of MacMinn and Primack (1996). 

The extragalactic 7-ray background, which can be accounted for as a superposition 
of spectra of unresolved blazars, and which we have predicted to be concave between 
0.03 and 10 GeV (Stecker & Salamon 1996), should steepen significantly above 20 GeV 
owing to our estimates of extragalactic 7-ray absorption at moderate to high redshifts. 
Both the predicted concavity and steepening may be too subtle to detect with present 
data from EGRET. However, next generation 7-ray telescopes which are presently being 
designed, such as GLAST, may be able to observe these features and thereby test the blazar 
background model. We also discuss the possible observability of dark matter lines in the 
multi-GeV region, given our predicted extragalactic background. 

We thank Yichuan Pei and Matthew Malkan for very helpful conversations. 
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Figure Captions 

Figure 1: Plot of f/-band mass-to-light ratios M/L versus stellar system age T for 
various values of metallicity Z, taken from Worthey (1994), Table 5A. The metallicity 
values are given on the right-hand side of the figure as y = \og{Z / Zq), where Zq is 
solar metallicity. Filled (open) squares correspond to the tabulated data for y = +0.5 
{y = +0.25), filled (open) circles for y = 0.0 {y = —0.22), filled (open) triangles for y = —0.5 
{y = —1.0), filled inverted triangles for y = —1.5, and open diamonds for y = —2.0. 

Figure 2: Plot of the S-band mass-to-light ratios M/L versus stellar system age T for 
various values of metallicity Z. See the caption of Figure 1 for details. 

Figure 3: Emissivity as a function of redshift, calculated using Eq. ^, for three 
wavelength values, A =0.28, 0.44, and 1.0 fim, for a Hubble constant value of ho = 0.5. 
Note that the emissivity scales as /ig (c/. Eq. 4 and Section 3). Solid line curves are for the 
case where the metallicity correction factor (£ from eq. |^) is used; dashed lines give the 
emissivity when this correction factor is not included. The data from the Canada-French 
Redshift Survey (Lilly et al. 1996) are also plotted. 

Figure 4: The intergalactic radiation energy density from stars and QSOs as a function 
of wavelength for redshifts z of 0, 1, 2, and 3, for a Hubble constant value of = 0.5 (the 
energy density scales as ho; see Eq. 7 and Section 3). These densities are calculated with 
the metallicity correction factor, £, included. 

Figure 5: The same as Figure 4, except without the metallicity correction factor. 

Figure 6: The opacity r of the universal soft photon background to 7-rays as a function 
of 7-ray energy and source redshift. These curves are calculated with the metallicity 
correction factor included in the expression for stellar emissivity. As discussed in the text, 
these results are independent of the value chosen for h^. 

Figure 7: The same as Figure 6, except without the metallicity correction factor. 
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Figure 8: The effect of intergalactic absorption by pair-production on the power-law 
spectra of four prominent grazars: 1633+382, 3C279, 3C273, and Mrk 421. The sohd 
(dashed) curves are calculated with (without) the metallicity correction factor. 

Figure 9: The extragalactic 7-ray background energy spectrum from unresolved 
grazars. The top and bottom sets of curves correspond to point-source sensitivities of 10~^ 
and 2 x 10~^ cm~^s~^, respectively, for 7-ray energies above 0.1 GeV, corresponding to the 
approximate point-source sensivities of the EGRET and GLAST detectors respectively. 
Because the FSRQ luminosity fuction that we employ scales as (Dunlop and Peacock, 
1990), our calculated EGRB spectrum scales as (see Eq. 10 in Stecker and Salamon, 
1996). 

Figure 10. The dot-dash polygon shows the allowed region of expected 7-ray photon 
counts calculated for a Bino {B) annihilation line as a function of M^. Present accelerator 
and cosmological constraints are indicated by the labels on the sides of the polygon (see 
text). In the figure labels, the letter "B" represents the Bino {B) and the letter "1" 
represents the slepton {€). As an illustration of how to read the figure, the arrow within the 
polygon indicates the line fiux upper limit for a Bino of mass 100 GeV. An exposure factor 
of 1 m^sr yr and an energy resolution of 10% are assumed. We also show the background 
fiuctuation count rate appropriate to these parameters for the lower set of fiux curves {i.e. 
with and without the metallicity correction) shown in Figure 9 (see text). 
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